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ABSTRACT 

A global lopsidedness in the distribution of the stars and gas is common in 
galaxies. It is believed to trace a non-equilibrium dynamical state caused by 
mergers, tidal interactions, asymmetric accretion of gas, or asymmetries related 
to the dark matter halo. We have used the Sloan Digital Sky Survey (SDSS) to 
undertake an investigation of lopsidedness in a sample of ~25,000 nearby galax- 
ies {z < 0.06). We use the m = 1 azimuthal Fourier mode between the 50% 
and 90% light radii as our measure of lopsidedness. The SDSS spectra are used 
to measure the properties of the stars, gas, and black hole in the central-most 
few-kpc-scale region. We show that there is a strong link between lopsidedness 
in the outer parts of the galactic disk and the youth of the stellar population 
in the central region. This link is independent of the other structural proper- 
ties of the galaxy. These results provide a robust statistical characterization of 
the connections between accretion/interactions/mergers and the resulting star 
formation. We also show that residuals in the galaxy mass-metallicity relation 
correlate with lopsidedness (at fixed mass, the more metal-poor galaxies are more 
lopsided). This suggests that the events causing lopsidedness and enhanced star 
formation deliver lower metallicity gas into the galaxy's central region. Finally, 
we find that there is a trend for the more powerful active galactic nuclei (the 
more rapidly growing black holes) to be hosted by more lopsided galaxies (at 
fixed galaxy mass, density, or concentration). However if we compare samples 
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matched to have both the same structures and central stellar populations, we 
then find no difference in lopsidedness between active and non-active galaxies. 
Indeed the correlation between the youth of the stellar population and the rate 
of black hole growth is stronger than the correlation between lopsidedness and 
either of these other two properties. This leads to the following picture. The 
presence of cold gas in the central region of a galaxy (irrespective of its origin) is 
essential for both star-formation and black hole growth. The delivery of cold gas 
is aided by the processes that produce lopsidedness. Other processes on scales 
smaller than we can probe with our data are required to transport the gas to the 
black hole. 

Subject headings: galaxies: structure, galaxies: interactions, galaxies: general 



Introduction 



In the standard AC DM universe, galaxies continue to grow significa ntly at the current 



epoch, accreting an ave rage of roughly 5% of their mass per Giga-year (e.g.. lKeres et al.ll2005 



Kauffmann et al.ll2006l ). Minor mergers, tidal interactions with close companions, and asym- 
metric accretion of cold gas can all perturb the underlying structure of the d ark matter halo, 
the stars, and the gas in galaxies for time-scales of-order a Giga-year (e.g. 



2007 



1997 : lLevine fc Sparke 



Mapelli et al. 



1998 



2008 



Cox et al. Il2008h . 



Zarits kv &: Rix 

Jog||l999l:lKornreich et allbooi iBournaud etal]i2005, : .Pi Matteo et al. 



A characteristic observational signature of these types of non-equilibrium situations is a 
global asymmetry (lopsidedness) in the distribution of the stars and/or gas. Indeed, galaxies 
commonly exhibit asymmetric structures. Asymm etry in the global H I 21cm emission- line 



profil e is present in half or more of disk galaxies (IRichter fc Sancisil 1 1994 : iMatthews et al. 
19981 ). Lopsid ed distributions of t he H I gas are revealed by spatial niaps of galaxies in 

The stellar 



the field (e.g., ISwaters et al.l 119991 ) and in groups (lAngiras et al. 



2006 



2007) 



mass distributions of disk galaxies are often lopsided as well. IZaritsky &: Rixl (119971 ) showed 
that ~30% of a sample of 60 field spiral galaxies had a significant m = 1 azimuthal Fourier 
component in the stellar mass. iRudnick fc Rixl (119981 ) showed that lopsidedness in the stellar 
mass distribution was common in early-type disk galaxies. Simulations of spiral galaxies have 
shown that the typical measured a mphtudes of lopsidedn ess cannot be explained by internal 
dynamical mechanisms alone (e.g. 
above are required. 



Bournaud et al.ll2005l ). The external processes described 



It has long been recognized that tidal interactions and mergers can act as triggers 
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for star-formation (IToomre fc Toomrd Il972l : iLarson fc Tinsleyl Il978l ). The resulting non- 
axisymmetric time-dependent gravitational potential can transport angular momentum out 
of the gas and also lead to shocks in which the gas is compres sed and kinetic energy is 



ultimately converted in t o radiation and l ost to the system (e.g., iMihos &: HernquistI Il996 



Di Matteo et al. 1 120071 : ICox et al. I l2008l ). These processes lead to the inflow of gas, an 



incre ase in gas density (either globally or locally), and hence an increase in the star formation 
rate (lKennicuttlll99^ 



Li et al.l (j2007al ) have investigated the link between interactions and star formation in 
the largest sample of galaxies studied to date. They found a strong correlation between the 
proximity of a near neighbor (closer than 100 kpc) and the average specific star formation 
rate. This result pertains to the relatively early stages of an eventual merger. Using lop- 
sidedness as a probe would allow us in principal to measure the history of the enhancement 
in star formation (in a n time-averaged sense) all the way through to the post-merger phase. 
Zaritsky &: Rixl (119971 ) exploited this idea, and found a correlation between lopsidedness 
and excess B-hernd luminosity emitted by a young stellar population. iRudnick et al.l (120001 ) 
showed that both recent (the past 1 Gyr) and ongoing star formation are correlated with 
lopsidedness, and estimated that as much as 10% of stellar mass in typical galaxies can be 
forr ned in such even ts. However, these papers investigated very small samples compared to 
the IlI et all J2007ah study. 



m 



An additi o nal li nk b etween lopsidedness and star formation is indicated by the results 



Keres et al.l (120051 ) and iBournaud et al.l (120051 ). The former authors use SPH simulations 



to conclude that the star formation history of galaxies is primarily regulated by the accretion 
rate of cold gas, rather than by merg ers. The accretion of t his cold gas will not occur in 
a spherically symmetric fashion, and iBournaud et al.l (12005! ) argue that this will lead to 
long-lived lopsidedness in disk galaxies. 

Large samples are crucial to the investigation of the c onnection between lop sidedness 
and the recent star formation history. In our recent paper ( Reichard et al.ll2008l . hereafter 
Paper I) we reported on our analysis of multi-color Sloan Digital Sky Survey (SDSS) im- 
ages of ~25000 low-redshift {z < 0.06) galaxies. We confirmed that Lopsidedness in the 
stellar mass distribution is indeed very common, and quantified its distribution as a func- 
tion of the principal structural properties of the galaxies. We found that galaxies of lower 
mass, lower dens it y, and lower concentration are systematically more lopsided. Similarly, 
Kauffmann et al. J2003al lb[) showed that galaxies with lower mass, density, and concentra- 
tion have younger stellar populations on-average. Taking these results together, the causal 
connections between star formation, lopsidedness, and other galaxy structural properties are 
ambiguous. Unraveling this web of mutual correlations requires the careful analysis of a 
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large and homogeneous sample. 



The strong corr elation between the s tellar mass of a galaxy and the metallicity of its 
interstellar medium (ITremonti et al.l 120041 ) provides a po werful constraint on the chemical 



evolution of galaxies and the intergalactic medium (e.g., iDalcantonI 120071 ) . If lopsidedness 
in a galaxy has been recently induced by either a minor merger with a gas-rich low mass 
companion galaxy or the accretion of cold intergalactic gas, then one signature would be a 
decrease in the metallicity of the interstellar medium in lopsided galaxies compared to other 
galaxies of the same mass. This can be tested with a large sample of galaxies. 

Interactions and mergers are also believed to be an important mechanism for fueling 
the formation and growth of a central supermassive black hole. This idea dates back at 
least to Toomre & Toomre (1972), and is one of the cornerstones of contemporary models 
that attempt to explain the co-evolution of galaxies and b l ack holes within the contex t 



of the hierarchical build-up of structure (e.g., iHopkins et al.l 120071 : iDi Matteo et al.l 120071 ). 
However, the observational verification of this idea remains insecure. At high-redshift the 
data are still too sparse to be conclusive, whi le at low - redshif t there have been many past 
studies that have led to contradictory results. iLi et al.l (l2007bl ) have undertaken an analysis 
of the largest low-redshift sample to date, based on 10^ galaxies in the SDSS. They find a 
strong link between close companions and star formation, but none between close companions 
and AGN. One possibility is that the black hole growth occurs only late in a merger, after 
the companion galaxy has been captured. This idea can be tested using lopsidedness and 
a similarly large sample since lopsidedness can be used to trace later stages of interaction 
than that seen in pair studies. 

In this study, we use the large data-set described in Paper I to try to understand how 
lopsidedness may be linked to star formation, metallicity, and the AGN phenomenon. We 
describe our specific galaxy sample in §2] and summarize our measurements of the structural 
parameters, star formation indicators, metallicity, and AGN indicators. Next, in §3l we 
compare the correlations between star formation and structure with an emphasis on sepa- 
rating out the dependence of star formation on lopsidedness versus other galaxy structural 
properties. Next, in §1] we explore whether the residuals in the mass- metallicity relation 
correlate with lopsidedness, and test whether this is a fundamental correlation, or only a 
secondary one induced through a mutual dependence of metallicity and lopsidedness on 
galaxy structure. Finally, in ^ we compare lopsidedness and the structural and stellar 
population properties for the host galaxies of AGN, to test for a possible link between in- 
teractions/mergers/accretion and the growth of black holes. We discuss the implications of 
our conclusions in ® 
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2. Data 



For the convenience of the reader, we briefly summarize the data, our sample selection, 
and our methodology here. The details are discussed in Paper I. The sample of gal axies was 
taken from the Sloan Digital Sky Survey (lYork et al.l 120001 : IStoughton et al.ll2002l ). a large 
survey of photometric and spectroscop ic data across tt sr. of the north ern sky. The sample 
is derived from SDSS Data Release 4 (jAdelman-McCarthy et al.ll2006l ). The survey's dedi- 
cated 2. 5 m telescope (jGunn et al.ll2006l ) at Apache Point Observatory uses a unique CCD 
camera (jGunn et al.lll998f) and drift-scanning to obtain w-, q-, r-, i-, and ^:-band photometry 
dPukugita et all llQod : Iflogg et all I2OO1I : llvezic et al.l booi ISmith et all I2OO2I : iTucker et al. 



20061 ). The pixel scale is O'.'396/px. Piber spectroscopy is obtained using 3" fibers and results 
in wavelength coverage between 3800—9200 A at a resolution R = \/6X = 1850 — 2200. 

Lopsidedness is defined here as the radially averaged amplitude of the m = 1 azimuthal 
Pourier mode between the radii in the SDSS images enclosing 50% and 90% of the galaxy 
light (i?5o and Rqq respectively). In Paper I we showed that galaxies on-average exhibit 
virtually identical distributions of r- and i-band lopsidedness (and only slightly larger values 
in the (7-band). We showed that weakness of the color dependence of lopsidedness implied 
that lopsidedness primarily traces the asymmetry in the underlying stellar mass distribution. 
Since any spatial variations in mass-to-light ratio will be smaller at longer wavelengths, we 
use the z-band lopsidedness A\ throughout this paper. 

The calculation of lopsidedness is subject to a handful of observational systematic er- 
rors. Insufficient spatial resolution, low signal-to-noise, and highly inclined galaxies in images 
lead to systematically incorrect measurements of Ai. Galaxies subjected to these system- 
atic effects were eliminated by requiring the following observational properties: z < 0.06, 
2i?5o/(PSP PWHM) > 4, total S/N > 30 in the region between R50 and Rgo, and b/a > 0.4. 
These restrictions yield a sample of 25155 galaxies that we use in our analysis below. These 
systematic effects prevent a reliable determination of lopsidedness in the central region of the 
galaxies (where the star formation properties are calculated). However, it has been shown 
that the lopsidedness i n mergers is very high in the central as well as the outer regions 
J Jog fc Mavbhatelbood l 



We make extensive use of ot her galaxy structural parameters derived from the SDSS 



data (see iKauffmann et al.ll2003a] ). The galaxy stellar mass (M^,) is derived from the SDSS 



z— band luminosity and a mass-to-light ratio derived from the spectra. The stellar surface 
mass density /i* is defined as the mean value within the galaxy half-light radius (/i* = 
M^:/27iRIq). We also use the galaxy stellar velocity dispersion a*, and the concentration 
parameter C = Rgo/R^o (the latter serving as a rough proxy for Hubble Type, increasing 
from ~ 2 for late type galaxies to ~3 for early types). 
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Most of the quantities we use to characterize the properties of the stars, ionized gas, 
and the AGN are based on the SDSS spectra. The SDSS fiber diameter is 3", which for 
our redshift-hmited sample corresponds to a typical projected diameter of about 3 kpc. Our 
measures of lopsidedness are made using the annulus between the radii enclosing 50% and 
90% of the light, a region lying entirely outside the SDSS fiber. Thus, it is important to 
emphasize that we are essentially relating the lopsidedness of the outer part of the galaxy to 
the central properties of the galaxy. 

We use a variety of me asurements to provide in formation about the star formation 
history in this central region (IKauffmann et al.ll2003al ). The 4000 A break, -D4000; indicates 
the current luminosity-weighted mean stellar age. Young stellar populations produce almost 
no metal absorption just shortward of 4000 A, giving -D4000 ~ 1; while older populations 
have strong metal line absorption, and -D4000 in creases to ~2 f or an age of ~ 10^° years. We 
use the narrow definition of 1^4000 as defined by iBruzuall (Il983l ): the ratio of the flux density 
Fi, in the ranges 4050—4250 and 3750—3950 A. Reddening affects this narrow definition less 
than other, wider definitions. 

The B.S absorption-line index HSa reaches maximum strength in a simple stellar pop- 
ulation with an age of 0.1-1 Gyr when A- s tars d ominate the continuum. We use the HSa 
index as defined by IWorthey fc Ottavianil (119971 ) where an equivalent width is calculated 
between two pseudo continuum bandpasses. The combination of 1^4000 and HSa allows us to 
recognize objects that have undergone bursts of star formation within the past Gyr (galaxies 
having abnormally strong HSa for their value of -D4ooo)- 



Recently, I Wild et al.l (120071 ) have applied the methodology of Principal Component 
Analysis to the near-UV spectral region using synthetic spectra derived from a large li- 
brary of star formation histories. The amplitudes of the first two principal components 
correspond closely to higher signal-to- noise measurements of -D4000 {PCi- mean age) and the 
excess in HSa for a given D4QQQ {PC2: a post-burst parameter). We use these parameters in 
our paper. 

We also use the specific SFR {SSFR or SFR/M^: the s t ar for mation rate per solar 
mass of stars in the galaxy) as computed by iBrinchmann et al.l (120041 ) based on the nebular 
emission-lines. Note that this is essentially an instantaneous measure of the central star 
formation rate (averaged over ~ 10^ years), and is hence complementary to the indicators 
above which have much longer time-constants. 



^This definition of post-starburst is made without regards to the emission-line properties (unlike some 
other common definitions). This definition allows the inclusion of both AGN and post-starbursts with some 
residual star-formation. 
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As we do not have direct measurements of the global stellar ages, we use the age-sensitive 
Petrosian g — i color from the SDSS images K-corrected to a redshift z = 0.1. This approach 
has the disadvantage that the color can be affected by dust extinction as well as by the mean 
stellar age. 

Type 2 (obscur ed) AGN are rec o gnized and characterized in the S PSS sample using 
the methodology in iKauffmann et al.l (l2003d ) and iHeckman et al.l (120041 ) . This is based on 
the detection of the signature of the narrow emission-lines that are excited by the AGN. 
The extinction-corrected luminosity of the [O III]A5007 emission-line is taken as an indica tor 
of the AGN bolometric luminosity and black hole accretion rate (see I Wild et al.l 120071 for 
a discussion of the extinctio n correction). The AGN hosts are early-type galaxies with a 
significant bulge component (IKauffmann et al.ll2003d ). and so the black hole mass (Mbh ) 
is estimated fror n a^, assuming the relatio n determined for bulges in iTremaine et al.l ( 12002| ) . 
As discussed by IKauffmann et al.l (l2003d ). for Type 2 AGN the SDSS images and spectra 
are dominated by starlight and so the properties of the host galaxies can be measured in 
the same way as for non-AGN. The exception is that because the AG N contaminate the 



nebula r emission-lines, the adopted SSFR in these cases is taken from iBrinchmann et al. 
( 120041 ) based on the mean relation between -D4000 and SSFR measured for non-AGN. 



The gas-phase metallicities for the star forming galaxies (non-AGN) are determined 
fro m fitting an extensive l ibrary of model galaxy spectra to the strong nebular emission-lines 



see 



Tremonti et al.ll2004j for details). 



3. Star Formation History and Lopsidedness 



Our plan for this section is to start with the simplest measures of the relationship be- 
tween the star formation histor y and lopsidedness. However, since both lopsidedn ess (Paper 
I) and star formation history (IKauffmann et al.l l2003bl : IBrinchmann et al.l |2004| ) correlate 
separately with other galaxy structural properties, we follow this with a more sophisticated 
multivariate approach. 



3.1. Simple Correlations with Star Formation History 

We begin with the Petrosian g — i color of the entire galaxy. In Fig. 1 there is a clear 
trend of increasing lopsidedness with bluer g — i color, with the median lopsidedness rising 
from 0.05 to 0.25 as g — i drops from 1.4 to 0.1. This correlation establishes that there is a 
strong connection between the lopsidedness in the outer region of the galaxy, and the star 
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formation history of the galaxy on global scales (as shown by iRudnick et al.ll2000l ). In a 
future paper we will use GALEX plus SDSS photometry to relate lopsidedness to the global 
star- formation rate. 

In the present paper our emphasis is on the central properties of our galaxies, so we 
now turn to stellar age indicators measured within the fiber. Each of our measures has a 
somewhat different time dependence. We consider them first individually and then together. 



We plot the distribution of A\ against -D4000 in the upper right panel of Fig. [H iKauffmann et al 



(j2003bl ) showed that the local galaxy population is bimodal in the distribution of -D40005 di- 
viding into young, late-type galaxies (-D4000 < 1-5) and old early-type galaxies (-D4000 > 1-7). 
For the young galaxies, the lopsidedness increases strongly with decreasing age: the median 
rises from 0.1 (mildly lopsided) at -D4000 = 1-5 to over 0.3 (strongly lopsided) at -D4000 ~ 
1. Values this low for -D4000 correspond to extremely young stellar populations (strong star- 
bursts, younger than a few tens of Myr). Most older galaxies (1^4000 > 1.7) uniformly have 
rather small values of A\ < 0.1. 

In the lower left panel of Fig. [T], we plot the lopsidedness distribution as a function 
of H6a- There is a correlation between median A\ and H6a for nearly the entire range of 
HSa, with the median value for A\ rising from near 0.05 at H6a = —1 to 0.20 at H6a = 9. 
Because values of H6a higher than about 7 can only originate about 0.1 to 1 Gyr after 
a strong burst of star-formation, we can connect such post-starbursts to events that have 
produced significant perturbations to the stellar mass distribution. 

In the lower right panel of Fig. [T|, we show the distribution of lopsidedness as a function 
the specific star formation rate (SSFR). Above SSFR ~ 10^^^ year^^ (or below a mass- 
doubling time of M^/SFR = 10^ Gyr), there is a strong correlation between star formation 
and lopsidedness. Typical star-forming galaxies with SSFR ~ 10^^^ to ~ 10^^° year~^ are 
mildly lopsided (median value of about 0.1) while strong starbursts (SSFR ~ 10"^ year^^) 
are very lopsided {A\ > 0.2). The smallest SSFR that can be reliably measured from the 
spectra is about 10~^^ year^^. Such galaxies are almost uniformly symmetrical (median A\ ~ 
0.05). 
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Fig. 1. — Two-dimensional distributions of A\ and star formation properties: g — i color, 
4000 A break, H6 absorption index, and specific star formation rate. The distributions of 
lopsidedness as functions of these properties are overlaid in gray: median (solid line) and 
quartiles (dashed lines). Lopsidedness tends to increase with bluer color, younger stellar 
population, and higher SSFR. 
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To orient the reader, in Fig. [2] we show a montage of galaxies selected to cover the range 
in mass and lopsidedness spanned by our sample of star-forming galaxies. The lopsided 
galaxies shown {A\ < 0.23) are not extremely distorted like classical major mergers. 
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Fig. 2. — Twenty-five SDSS galaxies with increasing lopsidedness from bottom to top and 
increasing stellar mass from left to right. Each image is 30" x 30", or about 23 kpc x 23 kpc 
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As noted above, the two age indicators 1^4000 and H6a used together can identify those 
galaxies that have undergone a recent (< Gyr) starburst, and separa te them from the more 
common galaxies experiencing more-or-less c ontinuous star form ation (IKauffmann et al. I l2003ah 
The parameters PCi and PC2 — PCi from IWild et al.l (120071 ) provide similar information 
at higher signal-to-noise. In Fig. [3], we examine the joint dependence of lopsidedness on 
both -D4000 and H6a in the left panel and on PCi and PC2 — PCi in the right. These 
panels both show that while the main correlation with lopsidedness is between the overall 
luminosity- weighted mean age (-D4000 or PCi), there is also a weaker trend for bursty galaxies 
at fixed mean age to be more lopsided. This trend can be seen in both plots, as the regions 
color-coded by median lopsidedness slant from lower left to upper right (implying increasing 
lopsidedness with the increasing strength of the post-burst signature for a fixed mean age). 
This result is confirmed in Table [T] where we list the correlation coefficient of lopsidedness 
with PC2, and also the partial correlation coefficient for the residual dependence of lopsid- 
edness on HSa after removing the mutual dependence on -D4ooo- The sample size is sufficient 
to give these correlations a large statistical significance. 

The highest mean values of lopsidedness {A\ > 0.2) are found along the far left edge 
of these panels. This region lies along the l ocus of a starburs t evolving to a post-starburst 
descendant over a period of about 100 Myr (IWild et al.l 120071 ). Galaxies in the region occu- 
pied by older post-starbursts (100 Myr to 1 Gyr) are significantly less lopsided {A\ ~ 0.1 to 
0.2). This is reasonable, given that the characteristic dynamical time over which the galaxy 
should relax after a m erger or tidal perturbat ion is a few hundred Myr. Indeed, the minor 
merger calculations of iBournaud et al.l (120051 ) show that values for lopsidedness as high as 
~0.2 last for only a few hundred Myr. 



These results are qualitatively consistent with those in IWild et al.l (120071 ) who inves- 
tigated a sample of relatively massive, bulge-dominated SDSS galaxies (required to have 
fi^ > 10^'^MQ/kpc^ and a* > 70 km s~^). They also showed that the largest average val- 
ues of lopsidedness are found in galaxies lying along the locus of starbursts/post-starbursts 
younger than about 100 Myr. However, the average value of lopsidedness they found for ob- 
jects lying along this locus was only about 60% as large as what we find in the same region 
of Figure [3] for our sample. Since we make no cut on fi^, or M^,, this difference is presumably 
related to the strong inverse correlations we found in Paper I between A\ and both and 
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Fig. 3. — Median lopsidedness as a function of the star formation indicators and principal 
spectral components. Starbust and post-starburst galaxies are typically the most lopsided, 
while galaxies with little star formation are typically the most symmetric. 
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3.2. Star-Formation History: Lopsidedness vs. Other Structural Properties 

We have shown above that there is a clear correlation between lopsidedness and all 
our indicators of the star formation history of galaxies. However, in order to establish a 
true causal connection it is necessary to demonstrate that this correlation is not simply 
induced by the known strong correlations of the basic galaxy structural parameters (m ass, 



density, and concentration) with both star formation history (IKauffmann et al.ll2003bl ) and 
with lopsidedness (Paper 1). 

We begin our assessment of the interdependence between lopsidedness, star formation 
history, and galaxy structure in Fig. HI In the three left panels we show the joint depen- 
dence of the median value of the age parameter -D4000 on lopsidedness and each of the three 
structural parameters. To get a complementary view, the three panels in the right column 
show the joint dependence of the median A\ as a function of 1^4000 and the three structural 
parameters. The analogous relations using the post-burst diagnostic PC2 — PCi and SSFR 
are shown in Fig. O and Fig. O These figures refer to the star formation history in the central 
region of the galaxy (inside the fiber). In Figure [7] we show the same set of plots, but now 
using the Petrosian g — i color as an indicator of the global star formation history of the 
galaxy. Our result is similar to Woods & Geller (2007) who find that the star formation 
becomes more centrally concentrated in interacting galaxies as the pair separation decreases 
and the specific star formation rate increases. 

In all the plots the overall trends are broadly similar: the left-hand panels show that 
galaxies with larger values of lopsidedness and smaller values of mass, density, and concen- 
tration have younger stellar populationsEI In each case, for a fixed value of a structural 
property (mass, density, or concentration), the stellar population becomes younger as the 
lopsidedness increases. Pictorially, we see that the bands color-coding the stellar age have 
diagonal or horizontal orientations. Similarly, the right hand panels in all three figures show 
that galaxies with younger stellar populations and lower masses, densities, and concentra- 
tions are more lopsided. In each case, for a fixed value of a structural property (mass, 
density, or concentration), the galaxies become more lopsided as the stellar population be- 
comes younger: again, we see that the bands color-coding the lopsidedness have a diagonal 
or horizontal orientation. 

We conclude from this that there is a correlation between lopsidedness and the star 



^In the left hand panels of Fig. [6] we are unable to calculate a meaningful value for the median SSFR for 
galaxies lying in the region of low lopsidedness and high mass, density, or concentration. This is because the 
SSFR in most galaxies is too small to be reliably measured. We leave this region blank, but the low median 
SSFRs there (< 10^^^ year~^) are consistent with what is seen in the figures using the other age indicators 
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formation history that is independent of any mutual correlation on the other galaxy structural 
parameters. To quantify this result, we turn to a partial correlation analysis to reveal 
which correlations are the strongest and perhaps the most fundamental. Table [T] shows the 
correlation coefficients of various combinations of structural and star formation properties. 
Values of logj^Q SFR/M^ below —12 have been set to —12 (see above). The results in this table 
show that there are separate and significant correlations between the structural properties 
themselves (independent of lopsidedness or star formation history), between lopsidedness 
and the other structural parameters (independent of star formation history), and between 
star formation history and lopsidedness (independent of the other structural parameters). 
These results are fully consistent with the visual representations in the figures. 

The lower left panels in both Figures [S] and [H] show that the relatively symmetric galaxies 
with A\ < 0.15 have higher values of PC2 — PCi and SSFR for less concentrated galaxies. 
However, the trend is reversed for the highly lopsided galaxies with A\ > 0.2. These lopsided 
galaxies typically show higher SSFR and larger PC2 — PCi at higher concentration. The 
most natural interpretation is that these are objects that are undergoing or have recently 
undergone central starbursts (leading to an increased concentration of light). This is con- 
sistent with the fact that an analogous behavior is not clearly seen in Figure [7] where we 
are investigating the global star formation history using the Petrosian g — i color. Here, the 
youngest (bluest) galaxies are lopsided but have low concentration. 
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Fig. 4. — Relationships between stellar age and structure shown in several slices of this 
parameter space. The color-coding of these two-dimensional histograms indicates the median 
-D4000 or A\. Galaxies with higher lopsidedness and smaller values of the other structural 
properties tend to have younger stellar populations. 
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Fig. 5. — Relationships between PC2 — PCi and structure shown in several slices of this 
parameter space. The color-coding of these two-dimensional histograms indicates the median 
PC2—PC1 or A\. Galaxies with higher lopsidedness and smaller values of the other structural 
properties tend to have greater values of PC2 — PCi and thus indicate recent star formation. 
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Fig. 6. — Relationships between specific SFR and structure shown in several slices of this 
parameter space. The color-coding of these two-dimensional histograms indicates the median 
SSFR or A\. Galaxies with higher lopsidedness and smaller values of the other structural 
properties tend to have higher SSFR, though the relationship with concentration is more 
complex. 
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Fig. 7. — Relationships between g — i color and structure shown in several slices of this 
parameter space. The color-coding of these two-dimensional histograms indicates the median 
g — i OT A\. Galaxies with higher lopsidedness and smaller values of the other structural 
properties tend to have bluer colors. 
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4. Metallicity and Lopsidedness 

The results above suggest that the processes that produce lopsidedness lead to an inflow 
of gas into the central region where it can fuel star formation. Any plausible source for this 
gas (capture of a small galaxy, gas transported from the outer disk to the inner region, the 
accretion of intergalactic gas) will have a lower metallicity than the pre-existing gas in the 
central region. Thus - provided that the inflow is rapid enough that the gas is not fully 
enriched by associated star formation before reaching the center - this inflow should have 
a signature in the chemical abundances in the interstellar medium. This has been seen by 
Kewley et al. (2006) and by Ellison et al. (2008a) for interacting galaxy pairs. Such a 
signature would also be consistent with Ellison et al. (2008b) who found that residuals in 
the mass-metallicity relation correlate inversely with the specific star formation rate. 

To test this idea, in Figure HI we show the mass-metallicity relation for our sample, 
color-coded by the median value for the lopsidedness in each cell. The strongest trend in 
this figure is simply the increase in lopsidedness with decreasing galaxy mass (Paper I). 
However, there is a weaker residual trend for the lopsidedness at a given value of stellar mass 
to increase with decreasing metallicity. This can also be seen in Table [2] where we list the 
relevant partial correlation coefficients. 



Tremonti et al.l (120041 ) showed that residuals in the mass-metallicity relation correlate 
with the surface mass density (with lower metallicity at fixed mass for less dense galaxies). 
We saw in Paper I that there is a strong inverse correlation between surface mass density 
and lopsidedness. Could the apparent correlation between lopsidedness and metallicity in 
Figure H] be induced through mutual correlations with surface density? The lower panel 
in Figure H] shows that at fixed surface mass density galaxies that have lower metallicity 
are more lopsided. This is confirmed by the partial correlation coefficients in Table [2j We 
therefore conclude that there is real correlation between lopsidedness and metallicity in the 
sense expected for an inflow of lower metallicity gas. 
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Fig. 8. — Relationships between metallicity and structure shown in two shces of this pa- 
rameter space. The color-coding of these two-dimensional histograms indicates the median 
A\. While the main correlation between the three parameters is between lopsidedness and 
either other structural parameter, there is a residual correlation between lopsidedness and 
metallicity at fixed mass. 
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We finish this section by quantifying the metalhcity deficit between the most and least 
lopsided galaxies of our sample. We use two subsamples taken from our main sample: low- 
lopsidedness galaxies {A\ < 0.08) and high-lopsidedness galaxies {A\ > 0.20). We compare 
the mass-metallicity relations of the two subsamples in Fig. |H The median metallicity of 
the less lopsided galaxies is 0.05—0.15 dex greater than in the more lopsided galaxies at the 
same mass. The metallicity differen ce is greater a t low m ass than at high mass. This deficit 
is consistent with the offset found in lEUison et al.l (120081 ) for galaxies with close companions. 
The relatively small amplitude of the effect would rule out an extreme picture in which very 
metal poor gas falls into the center and dominates the interstellar medium there. 




10.5 



Fig. 9. — Mass-metallicity relation for the most {A\ > 0.2, dashed lines) and least lopsided 
{A\ < 0.08, solid lines) galaxies. The 25th, 50th, and 75th percentile lines are shown for each 
distribution. The most lopsided galaxies have a metallicity deficit of 0.05-0.15 dex compared 
to the least lopsided galaxies at the same mass. The deficit is greater at low mass than at 
high mass. 
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5. Nuclear Activity and Lopsidedness 



We start with the simplest measures of the relationship betwee n AGN and lopsidedness. 
However, since both lopsidedness (Paper I) and AGN properties ( iKauffmann et al.ll2003c : 
Heckman et al.ll2004l : iKewley et al.ll2006bl ) correlate separately with other galaxy structural 
properties, we follow this with a multivariate approach to isolate the specific role of lopsid- 
edness. We also know that AGN with hig her luminosity are preferentially located in galaxies 
with younger central stellar populations (IKauffmann et al.ll2003d . 120071 ). and we have seen 
above that galaxies show a strong link between lopsidedness and a young central stellar 
population. We will attempt to sort out this complex web of correlations using an approach 
based on comparing AGN hosts to "twin" galaxies without AGN that have been matched 
to the AGN hosts in their structure and central stellar population. 

We begin with Figure [10] in wh ich we plot the d i stribut ion of lopsidedness in the di- 



agnostic line ratio dia gram used by IKauffmann et al.l (l2003c) to select Type 2 AGN (the 



BPT Diagra r n", aft er iBaldwin et al.lll98ll ). As discussed in 



Kauffmann et al.l (l2003d ) and 



Kewley et al.l (l2006bl ). the distribution running from the upper left to the lower right is the 
locus of star forming galaxies (with low metallicity at the top left and high metallicity at 
the lower right). The spur running to the upper right off the location of the metal-rich star 
forming galaxies is the locus of AGN. It represents a "mixing line" with AGN dominated 
objects at the upper right and star-forming/AGN composite objects at the lower left close 
to where the AGN spur joins the locus of star forming galaxies. 

It is clear from this plot that star forming gal axies are systemat i cally rn ore lopsided than 
AGN host galaxies. This is not surprising, since IKauffmann et al.l (l2003cl ) show that AGN 
hosts are generically massive, dense galaxies (which we have shown in Paper I to typically 
have small values for A\). In the region occupied by the AGN the median value of A\ is less 
than 0.1, with a weak trend for lopsidedness to be larger in the composite objects than in 
the pure AGN. This trend is not surprising, given the connection between lopsidedness and 
star formation discussed above. 

The Type 2 AGN in our sample span a large range in luminosity, and it is quite pos- 
sible that the mechanism responsible for fueling the growth of the black hole might change 
systematically as a function of AGN luminosity. To examine this, in Figure [TT] we plot the 
distribution of A\ vs. the ratio of the extinction-corrected [O III] 5007 emission line luminos- 
ity to the black hole mass. This ratio wil l be roughly proportio nal to the black hole accretion 
rate relative to the Eddington limit (see lHeckman et al.l 120041 ). 



We see in this plot that there is a systematic increase in lopsidedness with increasing 
AGN luminosity (increasing black hole growth rate). Below values of L[0 111]/ Mbh ~ 
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1Q~'^Lq/ Mq (corresponding roughly to accretion rates less than 10~^ of the Eddington limit), 
the median value of A\ is only ~ 0.05 (lopsidedness is essentially undetectably small). The 
median lopsidedness then increases with increasing AGN luminosity, reaching a value of 
~0.12 at L[0 m]/MBH = IOLq/Mq (~ 10% of the Eddington limit). This median value 
corresponds to only mild lopsidedness and (even at such high AGN luminosities) only about 
10% of the AGN hosts are strongly lopsided {A\ > 0.2). Thus, few of the high luminosity 
AGN hosts appear to be undergoing major mergers, but many are mildly disturbed. To give 
the reader a visual impression of the structure of AGN hosts we show a montage of SDSS 
color images of low- and high-luminosity AGN in Figure | 



Figure [TT] by itself does not establish a direct conne ction between AGN lumin osity 



and lopsidedness. As discussed in iKauffmann et al.l (l2003d ) and iKewley et al.l (l2006bl ). the 



properties of AGN hosts vary systematically as a function of AGN luminosity, and we have 
seen that lopsidedness is in turn strongly linked to galaxy structural properties (Paper I). 
Thus, in Figure [13] we examine the interdependence between lopsidedness, galaxy structure, 
and AGN luminosity. In the left three panels we see that, at fixed galaxy mass, density, and 
concentration, the median AGN luminosity increases with increasing lopsidedness. Similarly, 
the right-hand panels show that at fixed values of galaxy mass, density, and concentration 
the median value of lopsidedness increases with increasing AGN luminosity. The trends are 
most clear in the panels involving mass, and least clear in those involving concentration. 
These results are confirmed by the partial correlation coefficients listed in Table [31 

We conclude that there is a link between AGN luminosity and lopsidedness that is 
independent of the other galaxy structural parameters. This is analogous to the results for 
star formation, but the trends for the AGN are significantly weaker (the range spanned by 
the variation in the median lopsidedness in the AGN hosts is much smaller than in the star 
forming galaxies). 

The age of the stellar population in AGN host galaxies de c rease s systematically as a 



function of increasing AGN luminosity IKauffmann et al.l (12003d . l2007l ). We have also seen 



above that there is a very strong inverse correlation between lopsidedness and stellar age. 
Is the trend for more powerful AGN to be hosted by galaxies that are more lopsided simply 
induced by these other correlations? One way to visualize this is to plot the joint dependence 
of AGN luminosity on both stellar age (-D4000) and lopsidedness. This is shown in Figure [TH 
We see clearly that the primary correlation is that between the age of the stellar population 
and AGN luminosity with little dependence on lopsidedness (the color-coded bands of AGN 
luminosity are nearly horizontal). This is further quantified by the correlation coefficients in 
Table [3 



A more direct way to ask the question is as follows: for fixed values of both galaxy 
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structural and stellar age parameters, is there a remaining trend for higher luminosity AGN 
to be hosted by more lopsided galaxies? The large number of parametric dimensions involved 
in such a test leads us to adopt a "twinning" strategy. That is, for each AGN with a given 
value of L[0 111] /Mbh we will find within our sample another non-AGN galaxy with the 
same redshift, stellar mass, surface mass density, and stellar age (same value of -D4ooo)- The 
matching tolerances in redshift, mass, surface mass density, and 1^4000 were 0.005, 0.05 dex, 
0.05 dex, and 0.10, respectively. 

In Figure [151 we compare the distributions of lopsidedness as a function of the AGN 
luminosity for the AGN hosts and their non-AGN twin galaxies. In this plot, we have 
assigned each non-AGN twin the value for L[0 111]/ Mbh of its AGN twin. We have also 
plotted the distribution of the difference in lopsidedness (AGN host minus twin) as a function 
of AGN luminosity. We see from this figure that there is no systematic difference between 
the AGN hosts and their twins. We conclude that the trend for more powerful AGN to be 
hosted by galaxies that are more lopsided is due to the strong link of the age of the stellar 
population in the galaxy bulge to both the AGN luminosity and to the lopsidedness of the 
host galaxy. This can also be seen in Figure [12] where we show the very similar dependence 
between stellar age and lopsidedness for both the AGN hosts and their twins. These results 
are quantified by the partial correlation coefficients in Table [31 

This leads us to the following thoughts. It is clear that the processes that produce 
lopsidedness facilitate the delivery of gas into the central regions of galaxies. It is also clear 
that the presence of cold gas in the central region (however it arrives there) is required for 
the formation of stars there. The fact that the more rapidly growing black holes are strongly 
associated with a younger stellar population implies that the presence of cold gas in the 
central region is also required for rapid fueling of the black hole. The results above then 
suggest that the once the gas has arrived in the central kpc-scale region, some other process 
(not directly related to the cause of the lopsidedness) regulates the subsequent transfer of a 
small fraction of the gas inward by orders-of-magnitude in radius to the black hole accretion 
disk. 



It is interesting to compare our results to those of lLi et al.l (l2007bl ). who have examined a 
very similar SDSS galaxy sample for a link between the presence of close companion galaxies 
and the fueling of black holes. They found that close companions are associated with a 
significant boost in the star formation rate, but are not associated with an enhancement of 
nuclear activity (compared to galaxies of similar mass). Ellison et al. (2008a) reached similar 
conclusions with a different sample drawn from SDSS. These results contrast with our result 
showing that (at fixed galaxy mass, density, or concentration) higher AGN luminosities are 
preferentially found in galaxies that are more lopsided. 
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Li et al.l (j2007bl ) concluded from their results that the star formation induced by a close 
companion and the star formation associated with black hole accretion are distinct events. 
They argued that these events may be part of the same physical process, for example a 
merger, provided they are separated in time. They therefore speculated that accretion onto 
the black hole and its associated star formation occurs after the two interacting galaxies have 
merged. This scenario is consistent with our results, since lopsidedness will persist for several 
galaxy dynamical times after the merger is complete. In a future paper we will combine our 
results on close companions and lopsidedness to try to test this idea. 
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Fig. 10. — BPT diagram color-coded by median A\. The strong correlation with lopsidedness 
is that with star formation rather than that with AGN acti vity. The black line sepa rates 
the AGN from the star-forming galaxies and is adopted from iKauffmann et al.l (j2003d ). 
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Fig. 11. — The distribution of A\ and L[0 IH]/Mbh for all galaxies in the sample. The 
5th, 25th, 50th, 75th, and 95th percentiles of A\ are overplotted in gray. Typical lopsided- 
ness increases with L[0 HII/Mbh- Median A\ reaches a moderate value {A\ = 0.1) when 
L[0 111]/Mbh is high at 10 Lq/Mq. 
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Fig. 12. — Twenty-five SDSS galaxies with increasing lopsidedness from bottom to top and 
increasing L[0 Vl[\IMbh from left to right. Each image is 30"x 30", or about 23 kpc x 23 
kpc at z = 0.04, a typical redshift of the sample. 
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Fig. 13. — Relationships between L[0 IH]/Mbh and structure shown in several slices of 
this parameter space. The color-coding of these two-dimensional histograms indicates the 
median L[0 IU]/Mbh or A\. In each case, there is a correlation between lopsidedness and 
AGN luminosity independent of these structural parameters. 
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Fig. 14. — Relationship between L[0 IU]/Mbh, -D4000, and structure. The color-coding of 
these two-dimensional histograms indicates the median L[0 111]/ Mbh- The primary corre- 
lation is between stellar age and L[0 111]/ Mbh- There is no noticeable correlation between 
L[0 111]/ Mbh and lopsidedness in this data. 



- 31 - 




-3-2-1012 
log,o {L[OIII]/M,h) 



Fig. 15. — Top left and top right: Two-dimensional distributions of L[0 111]/ Mbh and 
lopsidedness of non-AGN and AGN galaxy pairs matched in redshift, mass, mass density, and 
stellar age (-D4ooo)- The distributions show no significant difference between the lopsidedness- 
AGN luminosity relationship. Bottom: Difference in lopsidedness of the AGN and non-AGN 
in matched galaxies pairs vs. the AGN luminosity. No preference is shown for the AGN 
or non-AGN to be more often more lopsided than the other for any value in our range of 
L[0 m]/MBH- 
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Fig. 16. — Relationship between lopsidedness and mean stellar age (-D4000) for the AGN and 
non-AGN in pairs of galaxies matched in redshift, mass, mass density, and mean stellar age. 
The star formation-lopsidedness connection is similar for AGN and non-AGN. 
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6. Summary 



We have studied a sample of approximately 25000 low-redsliift {z < 0.06) galaxies from 
the Sloan Digital Sky Survey (SDSS) Data Release 4 to investigate the links between large- 
scale asymmetries in the stellar mass distribution in galaxies (lopsidedness) , star formation, 
the metallicity of the interstellar medium, and the presence of active galactic nuclei (AGN). 
Lopsidedness has been defined as the radially averaged m = 1 azimuthal Fourier amplitude 
(Ai) measured between the radii enclosing 50% and 90% of the galaxy light in the SDSS 
images (i.e. in the outer part of the galaxy). We have previously shown that lopsidedness 
traces the distribution of the underlying stellar mass (Paper I). Lopsidedness is a signpost 
of a non-equilibrium global dynamical state, and can be induced by mergers, asymmetric 
accretion of cold gas, tidal interactions, or underlying asymmetries related to the dark matter 
halo. We have used spectra obtained through the SDSS 3" diameter fibers (typical projected 
diameter of ~ 3 kpc) to characterize the stars, gas, and AGN in these central regions. 

We have used the amplitude of the 4000 A break, the strength of the high-order Balmer 
absorption-lines, and the specific star formation rate derived from the nebular emission- 
line to characterize the stellar population in the central few-kpc-scale region. We found 
strong links between lopsidedness and recent /on-going central star-formation: galaxies with 
younger stellar populations are more lopsided and more lopsided galaxies have younger stellar 
populations. Starburst and post-starburst galaxies are the most lopsided on average. 

We have previously shown that galaxies with lower surface mass density and mass 
are more lopsided (Paper I) , and that galaxies with lower mass and lower d ensity have 
younger stellar populations (IKauffmann et al.l l2003bl : iBrinchmann et al.l l2004j ). Here, we 
have shown that there is still a strong correlation between lopsidedness and the age of the 
stellar population even after their mutual dependences on these other structural parameters 
have been removed. These results are consistent with other evidence that mergers and tidal 
interactions trigger central star-formation, but place these results on a firm statistical base. 
They are also consistent with the idea that star formation in galaxies toda y is regulated by 
the a ccretion of cold gas (IKeres et al.ll2005l ) which can excite lopsidedness (IBournaud et al. 
20051 ). They imply that the timescale for the transport of gas into the central region can not 
be significantly longer than the timescale over which lopsidedness persists in the outer disk. 
This is consistent wi th recent numerical siraulations of star fo rmation in minor mergers and 
galaxy interactions (IDi Matteo et al. 1120071 : ICox et al. II2008I ). 



We have also shown that at fixed stellar mass, more lopsided galaxies have system- 
atically lower gas-phase metallicities by about 0.1 dex on average. This suggests that the 
processes causing lopsidedness deliver lower metallicity gas into the central region, but the 
small amplitude of the effect rules out extreme events as being typical. 
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Finally, we found that there is a trend for the more rapidly growing black holes to be 
hosted by galaxies with higher average lopsidedness. This is true when comparing galaxies 
at fixed galaxy mass, surface mass density, and concentration. However, when the AGN 
hosts were compared to non-active galaxies that were also matched in the age of the stellar 
population in their central region, we found no excess lopsidedness in the AGN hosts. The 
strongest link is that between the youth of the stellar population and the growth rate of 
the black hole. The correlations of these two properties with lopsidedness are weaker. This 
suggests that the presence of cold gas in the central few kpc-scale region (which is facilitated 
through the processes that produce lopsidedness, and which leads to significant star forma- 
tion) is a necessary but not sufficient condition for subsequent fueling of the growth of the 
black hole. Other processes are subsequently required to deliver the gas from scales of a few 
kpc all the way to the black hole accretion disk and these would not be directly related to 
the process(es) that produced lopsidedness. 



Combining our results with the analysis by iLi et al.l (l2007bl ) and Ellison et al. (2008a) 
of the role of close companion galaxies in driving star formation and fueling black holes, 
suggests that the period of black hole growth may be preferentially associated with the end 
stages of a minor merger. This idea will be tested in a future paper. 
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Table 1. Partial Correlation Coefficients: Lopsidedness and Star Formation 
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Table 2. Partial Correlation Coefficients: Lopsidedness and Metallicity 
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Table 3. Partial Correlation Coefficients: Lopsidedness and AGN Activity 
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